Magnetic layers with periodic point 
perturbations 

p. Exner/'^ K. Nemcova^'^ 



1 Nuclear Physics Institute, Academy of Sciences, 25068 Rez 
near Prague, Czech Republic 

2 Doppler Institute, Czech Technical University, Bfehovd 7, 
11519 Prague, Czech Republic 

3 Institute of Theoretical Physics, FMP, Charles University, 
V Holesovickdch 2, 18000 Prague, Czech Republic 

exner @uj f . cas . cz , nemcova@uj f . cas . cz 



Abstract 

We study spectral properties of a spinless quantum particle confined 
to an infinite planar layer with hard walls, which interacts with a peri- 
odic lattice of point perturbations and a homogeneous magnetic field 
perpendicular to the layer. It is supposed that the lattice cell contains 
a finite number of impurities and the flux through the cell is rational. 
Using the Landau-Zak transformation, we convert the problem into 
investigation of the corresponding fiber operators which is performed 
by means of Krein's formula. This yields an explicit description of the 
spectral bands, which may be absolutely continuous or degenerate, 
depending on the parameters of the model. 

Keywords: Schrodinger operator with magnetic field, Dirichlet layer, 
periodic point potential, group of magnetic translations 

1 Introduction 

The object of the present study is a three-dimensional spinless quantum par- 
ticle interacting with a homogeneous magnetic field and a periodic point po- 
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tential. In addition, the particle is confined to a fiat layer of a constant width 
d with Dirichlet boundary conditions, the magnetic field being perpendicular 
to the layer. The lattice of the point interactions is two-dimensional, not 
necessarily planar, and both of its generating vectors are, of course, parallel 
to the boundary planes of the layer. 

The problem is clearly motivated by the need of modeling electrons in a 
semiconductor layer with 'impurities', either natural or artificially created. 
If they may be supposed to be well localized, one often describes them by 
point interactions. It is useful because in such a way we get solvable models 
which make it possible to derive the spectral and scattering properties of a 
given configuration in a relatively simple way, and at the same time they 
reproduce the basic features of an actual crystal-lattice layer with some alien 
atoms and a low-density electron gas. Various systems of this type have been 
investigated in the literature; we refer to our previous paper ^Sj and an earlier 
study of the two-dimensional analogue ^ for an extensive bibliography as 
well as for a discussion of the used approximations. 

The said paper was devoted to investigation of a flat hard-wall layer 
with a finite number of point interactions, with or without a magnetic field; 
we analyzed there the spectra, and in the former case also the scattering 
properties of such systems. The situation is more complicated when the 
number of point perturbation is infinite. The spectral content of such mod- 
els is very rich and difficult to treat in the general setting; this is why we 
restrict our attention here to the particular case specified above, in which 
the periodicity allows us to employ the Bloch-Floquet decomposition, with 
its specific features due to presence of a homogeneous magnetic field that 
were first pointed out by J. Zak in his classical papers [T^ . 

On the other hand, the present work represents a generalization of the 
results derived in [7] for a two-dimensional Schrodinger operator with a pe- 
riodic point potential and a homogeneous magnetic field perpendicular to 
the plane. Loosely speaking we 'add' the third dimension in the direction of 
the magnetic field and allow the positions of the point potentials to 'spread' 
transversally preserving the periodicity - this is what we meant by a two- 
dimensional lattice in the first paragraph. 

Our basic tool for spectral analysis is the Krein formula which expresses 
the resolvent perturbation due to the point interactions. It is important that 
in absence of the impurities the Hamiltonian of the system allows a separation 
of the plane variables from the transverse one. As a consequence, the free 
resolvent kernel of the fiber operator (as well as all quantities derived from 
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it such as Bloch eigenf unctions, etc.) can be written by means of explicitly 
given series. In this sense our model is solvable, the only difference from 
those in the full space is that there such quantities can be written in terms 
of elementary or special functions. 

In accordance with [7j we adopt two more simplifying assumptions: we 
suppose that the number of point interactions in an elementary cell of the 
lattice is finite and the flux through the cell is rational. Especially the second 
one is important because one naturally expects in analogy with the two- 
dimensional case that the spectral character can be substantially different 
for an irrational flux; it is an independent problem of its own interest which 
we are not going to discuss here. On the other hand, following the setting of 
the paper [7j we study general point perturbations including the situations 
when the particle can "jump" between different impurity sites within a cell. 
The reason is that mathematically it does not mean a lot of extra work; 
physically interesting case is naturally that of point interactions defined by 
local boundary conditions, i.e. with the parameter matrix A appearing in 
the relation (j4.H) being diagonal. 

Let us review briefly the contents of the paper. We start in the following 
section with analyzing the free Hamiltonian. Next, in Sectional we formulate 
the Bloch-Floquet theory for the present case; the main result is so-called 
Landau-Zak transformation which makes it possible to reduce the spectral 
problem to investigation of suitable fiber operators. This is done by means of 
Krein's formula in Section^] The last three sections are devoted to successive 
analysis of the spectra at three levels of complexity: (i) an integer flux (in the 
units of magnetic-fiux quanta) through the elementary cell and a single point 
interaction in the cell, (ii) an integer flux and a finite number of impurities 
in the elementary cell, (Hi) and finally, the case of an arbitrary rational flux. 

2 Free Hamiltonian 

Consider an infinite layer of a fixed width d, i.e. S = x [0,d]. The 
coordinates we are going to use are x = {x,X3), where x = {xi,X2) G 
and X3 G [0,d]. The layer is placed into a homogeneous magnetic field of 
intensity B = {0,0, B). We will use the circular gauge, A = X2, Xi, 0), 
and rational units, h = c = e = 2m = 1. A quantum particle confined to the 
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layer is described by the free magnetic Schrodinger operator in L^(S), 

Ho = i-tW-A)^ (2.1) 

with Dirichlet boundary conditions 

^{x,0) = 4j{x,d) = 0, xeR'^. (2.2) 

In the absence of an additional interaction the variables separate, i.e. the 
operator Hq can be decomposed into transverse modes, 
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where I is the unit operator in L'^{0,d). The first two terms at the r.h.s., 
in the following denoted as h, describe a two-dimensional particle in the 
perpendicular homogeneous field in the circular gauge. The third term rep- 
resents the energy of the nth transverse mode; the corresponding normalized 

eigenfunction will be denoted as Xn ■ Xni^^) = sin (^^^). The resolvent 

kernel of the operator h is well known [3j to be 

{h — z)"^ {x,x') = — exp ( —i— X A x' — ^—^ \x — x'\'^ 
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where U is the irregular confluent hypergeometric function, T is the gamma 
function and x A X . — X'^OC'2 OC'^'JC'^' For the sake of brevity we denote this 
kernel as Gq^(x, x'; z). The decomposition ()2.3p then yields an explicit form 
for the resolvent kernel of the operator (j2.H) . namely 

G'o(x, x';z) = {Hq- zY^ {x,X3]x\x'^) 
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where kn a/ -2 — (nn/d)'^. Applying Fourier transformation to the trans- 
verse part, M: L'^{[0,d])^ e'^(N) defined by 

{Mf){n)^ f dyXn{y)f{y), neN, (2.5) 
one can rewrite the free resolvent kernel as a matrix with the elements 

Gq{x, n; x', n'; z) = Sn„'Gl^ (^x, x'; z - (^^^ ^ . (2.6) 

Consequently, the spectrum of the free Hamiltonian consists of Landau levels 
shifted by the energies of the transverse modes, 

a{Ho)=aess{Ho) = i^\B\{2m + l) + (^^y : m,n-leNo|. (2.7) 

The energy corresponding to the Ith Landau level and nth transverse mode 
will be denoted as e{l,n) := \B\{21 + 1) + {nn/dy. If the \B\ and (n/d)'^ 
are rationally related it may happen that there exist more than one pair of 
numbers {li^rii) giving the same value, e{li,ni) = Zq. In such a case we will 
denote this set of pairs as J{zo) and the number of these pairs as | J{zo) \ ■ 



3 Magnetic translations 

Next we consider a lattice T = A + K, where A is a lattice spanned by two 
independent vectors in M^, which can be without loss of generality supposed 
to be a = (ai, 0, 0) and b = (6i, 62, 0) with 62 7^ 0, and i^' is a set of \K\ points 

i = 1, . . . , n, from the elementary cell {sa + tb : s, t G [0, 1)} x (0, d) of 
the (interior of the) layer. 

As usual with periodic systems, the first thing to do is to find the ap- 
propriate representation of the translation group. In presence of a magnetic 
field the argument shift must be composed with multiplication by a suitable 
phase factor. Specifically, the translation by a vector v = (w, 0) acts in the 
state space -L^(E) of our system as follows, 

f{x)>—>-exp{—ni^xAv)f{x—v), (3.1) 

where ^ is the number of the flux quanta of the fleld B through the unit 
area of the plane M^. Recall that the quantum of magnetic flux is given 
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by 27r^c/|e|, thus with the chosen system of units we have ^ = B/{2tt). 
The phase factor at the r.h.s. can be rewritten as exp {^{B x v).S^. Since 

B = (0, 0, B) and we consider translations in the plane, i.e. by a vector 
V = (fi,f2,0), it is obvious that the third coordinate plays a marginal role 
only in magnetic translations. This allows us to follow closely the procedure 
used in |7j to analyze spectral properties of a particle confined to a plane, 
i.e. with the third components of the vectors v and x absent. 

Let us briefly summarize needed notions and facts. First we define the 
group of discrete magnetic translation over the lattice A, 

W{^, A) = {(A, C) : A G A, C = exp(7rir/n), n G Z} , (3.2) 

where rj = 0162^ = aAbC, is the number of flux quanta of the field B through 
the elementary cell. In the coordinates relative to the basis a and b the 
multiplication in W{C,, A) has the form 

(A,C)(A',C') = (A + A',CC'exp(7rzr^(A,A;- A,A1))) . (3.3) 

Notice that groups W{^i,Ai) corresponding to different values of and dif- 
ferent lattices Aj but having the same value of t] are isomorphic, hence we 
will denote the group simply by Wrj. Next we define the representation T of 
the group Wrj in the space L'^iT,) as follows 

(r(A,C)/)(x) =Cexp(-7r^exAA)/(f-A). (3.4) 

Furthermore, replacing s G S in the last formula by vectors 7 G F we get a 
representation D of Wrj in the space i'^(F). 

If the flux ?7 is a rational number, then any unitary representation of the 
group Wrj can be uniquely decomposed into orthogonal sum of irreducible 
representations. Any 'physical' irreducible representations can be parame- 
terized by a point p = (pi,P2) from the torus = [0, 1) x [0, 1) - see [TT]. If 
7] is an integer, then the group is Abelian and the irreducible representations 
are one-dimensional, their characters being given by 

X((A,C);P) = Cexp ( - 2m{XaPi + Xm + {N/2)K\,)) . (3.5) 

For a general rational flux, t] = N/M, the irreducible representations of 
are generally M-dimensional. In particular, the generators (a, 1) and (6, 1) 
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are in this case given by 



A((a, l)-p) = diag [e-2-*Pi,e-2"^(Pi+''), . . . ^ e-2"*(P^+(*^~^)'')] 

hi-i 

g-27rip2 g 



A((6,l);p) 



(3.6) 



where zeros in the second hne are (M— l)-dimensional column and row vectors, 
respectively. In order to obtain nonequivalent representations, of course, we 
must restrict here p to the torus = [0, M~^) x [0, 1). 

The decomposition into irreducible representations is accomplished by the 
Landau-Zak transformation Crj. Recall that the magnetic translations from 
Wr] does not affect the third coordinate, or in the transverse-modes repre- 
sentation, they do not affect the quantum number n. By a straightforward 
modification of the two-dimensional formulae [7j, we define thus C^j by 



(>Cr,/)(p, j, = A^"5 exp hnim ^ ^ j 

irtn — rv~i ^ ' 



(3.7) 



m=— oo 

X / cfxf{x,n)%jjQ{x;pi + rij + m,l) 

M—lN-1 oo oo pl/M pi 



iq'f)ix,n)=N''^J2J2T.Yl / / dp2fip,J,kJ,n) 



j=0 k=0 1=0 m=—oo 

—27Tim — — — j iIjq{x;pi + r^j + m, /) , 

where '?/'o(a;; 5', /) with g G M and / = 0, 1, ... , are generalized eigenfunction 
of the operator h (the 'planar' part of the Hamiltonian Hq) associated with 
the lattice A, 

«-')-(^^^'-«)"'-("^^0-H(ef-')) 

X exp \B\(^X2 + ^ ^ Hi(^\B\^ (^X2 + ^ q 

with Hi being the /th-order Hermite polynomial. As in the two-dimensional 
case, we get the following claims. 



7 



Theorem 3.1 Assume that the flux is a rational number, rj = N/M. Then 
the Landau-Zak transformation j.^. ?| ) decomposes the representation T of a 
magnetic translation {v, () G Wj^ into a direct integral of multiples of irre- 
ducible representations A(-), in other words 

Cr,T{vX)jC-' = / d^pA{vX;p)(»Ic^®hHN)®h^N). (3.8) 

In particular, if the flux rj is integer, then the irreducible representations are 
one-dimensional and their characters are given by the expression l[y.5]) . 

The second important feature of the Landau-Zak transformation is that it 
diagonahzes the free Hamihonian in the space L^(T^), as weU as any other 
operator which exhibits the magnetic-translation symmetry: 

Theorem 3.2 Assume that a self-adjoint operator H acting on i^^(S) is 
invariant w.r.t. the group Wr,; then it can be decomposed into a direct integral 

H= d^pH{p), (3.9) 

where H = [Crj ® M.)H{Cri ® M)^^ and the fiber operator H{p) acts on the 
space C*^ (g) ® » 



4 The perturbed Hamiltonian 

Now let us consider perturbation of the free Hamiltonian Hq by point po- 
tentials placed at the points of the lattice F. To define the perturbed 
operator H^, where A are the coupling-parameter matrix introduced by 
the relation ()4.1|) below, we employ the standard technique based on the- 
ory of self-adjoint extensions of symmetric operators. First, we pass to 
the symmetric operator which is the restriction of Hq to the domain 
D{Sr) := {/ G D{Ho) : /(f) = 0, 7 G F}, it is well defined in view of the 
usual Sobolev embedding. The sought Hamiltonian is then an appropri- 
ate self-adjoint extension of the operator Sr- There is a family of such exten- 
sions which can be parametrized in different ways, the most traditional one is 
the von Neumann method using unitary maps between deficiency subspaces. 
For our purposes it is more suitable to employ the mentioned self-adjoint 
operator A acting in i'^(F); its properties will be specified later. 
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One of possible characterizations of the extensions starts from the (gen- 
erahzed) boundary values. In this approach, one determines the domain 
D{Ha) as the set of all functions having prescribed behavior in the vicinity 
of the point-potential sites 7 G F, namely 

\x 71 

where the vectors Lj{ip) := {Lj{ip; 7) : 7 G F } , j = 0, 1 , are related by 

L.iij) = -AttALoW . (4.1) 

This is a generalization of the usual point-interaction definition - see |2] - 
the latter corresponds to a diagonal A and represents the physically most 
interesting situation, at least from the viewpoint of modeling semiconductor 
layers with impurities, as we pointed out in the introduction. 

Spectral properties of the Hamiltonian can be found by means of its 
resolvent. Krein's formula gives us its kernel, i.e. the Green function of Ha 
for a fixed operator A, 

Gix,x';z)=Go{x,x';z)- [Q{z) + A]-\i,i')Go{x,f, z)Go{f , x'; z) . 

7,7'Gr 

(4.2) 

with Q{z) given by the relation (j4.4j) below. If we use the transverse-mode 
representation of the Hamiltonian in the space L^(R^) ® £^(N) we get 



G{x,n; x' ,n; z) = 




Since the Krein's formula ()4.2j) was originally meant for a finite number of 
point potentials and the set F is infinite, we have to be a bit more precise. 
One possible way how to define the Hamiltonian Ha is via strong resolvent 
limit of the family of "restricted" operators ^^^(f ) which are point- interaction 
Hamiltonians referring to finite subsets F C F, an example for diagonal A 
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and 5 = can be found in [21 III. 1.1]. But the limit is not necessary; the 
existence of the operator Ha and the generahzed Krein formula have been 
already proven for a larger class of operators, - see, e.g., or 
The matrix Q in ()4.2p is given by 



Q(7,7';^) 



Co (7,7'; z) 



7 = 7' 



(4.4) 



where Qo(73; z) is the regularized Green function (stripped off the pole sin- 
gularity) which is defined by lim|^_^|^o(G'o(a;, 7; z) — 4^|x — 7|~^)- We know 
from [S] how it looks like. 
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In case when two point potentials are arranged vertically, i.e. 7 = 7' and 
73 7^ 73 5 the corresponding element of Q is well defined but the expression 
(12. 4j) makes no sense and has to be recast into the form 
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For the sake of brevity it is useful to rewrite the Krein's formula in a compact 
form, 

Ra{z) = Ro{z) -T,[Q{z) + A]-'t:, 
where is a map from f{T) to L^{R'^) (g) f(N) defined by 



(4.6) 



(r,0) ix,n):=J2 U 7; ^ - (^) 



Xn(73 



(4.7) 
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The operator Hq is invariant under magnetic translations from the group Wjj. 
One can easily check that the free resolvent kernel satisfies the relation 

Go{x — \,x — A; z) = exp(7ri^ (x — x') A \)Go{x, x; z). (4.8) 

Our aim is to study the situation when the operator is also PV^-invariant, 
which is the case when the operator A is assumed to satisfy the same con- 
dition. Notice that this is trivially satisfied if A is a diagonal matrix. We 
will assume only that A is a self-adjoint operator invariant w.r.t. VF^, and 
furthermore, that it is bounded and there exist two positive constants Ci and 
C2 such that 

1-4(7, 7')| < Ciexp(-C2|7 - 7'|) for all 7, 7' G T . (4.9) 

The last condition means a restriction on the non-locality we have allowed 
mathematically: it means that the probability of particle hoping between 
two points of the lattice F decays exponentially with their distance. 

A similar estimate is valid for the operator Q - for any given z ^ cr{Ho) 
there exist two positive constants C3 and C4 such that 

|Q(7,7';^)I <C3exp(-C4|7-7l) forall7,7'er, (4.10) 

as it follows from the definition ()4.4|) of the function Q and the free resolvent 
kernel ()2.4|) . The infinite sum contained in the second formula converges be- 
cause the term r{u)U{u, 1; s) can be written for large positive u as 2KQ{2y/us) 
(see [H 13.3.3]) and the Macdonald function Kq decays exponentially for large 
argument. Using the asymptotics of the function U{u, 1; s) we conclude that 
the sum grows with I7I at most as a polynomial; hence the exponential term 
exp(— I7P) is sufficient to yield the estimate. 

If z approaches a point from the spectrum of the free Hamiltonian Hq the 
elements of Q may diverge - cf. ^ - in other words the functions (5(7,7'; ■) 
are in general meromorphic. We summarize the above discussion in the 
following theorem: 

Theorem 4.1 Suppose that an operator A acting on £^(r) is self-adjoint 
and W^-invariant. Then there is exactly one self-adjoint extension Ha of the 
operator Sr with Green function given by 

G{x,x'-z) =Go{x,x'-z) - J2 [Q{z)+A]-\^,f)Go{x,f,z)Go{f,x'-z), 

7,7'er 
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where the operator Q{z) is defined by the relations ^■4\ ) o-nd \4-^ - As an 
operator in d-'^iT), i.e. an infinite matrix, (5(7,7'; 2) satisfies the estimate 
for some C3, C4 > and it is Wrj-invariant. As a function of z, Q 
is meromorphic and all its poles belong to the spectrum a{Ho). Finally, the 
operator Ha is also Wjj-invariant. 

5 The case of integral flux and a monoatomic 
crystal 

We begin with the simplest case assuming that the flux r] through the elemen- 
tary cell is integral, f] = N > 1, and that this cell contains only one potential 
placed at K = (0, 0, K3); in other words in this section we have F = A + {k}. 

We already know from Thm. 13.11 that to diagonalize the representation 
T acting in i^^(S) we have to employ the Landau-Zak transformation. Since 
M = 1 in the present case we drop the parameter j. Diagonalization of the 
representation D acting in £^(r) is achieved by the Fourier transformation 

■■ i\A-) ^ L\T^) which is defined by 

(^,(/))(p) = ^0(A + /t)eA(p) (5.1) 
XeA 

with the basis e\{p) := exp(— 27rz(AaPi + XbP2 + NXaXt/^)). It reduces the 
action of the magnetic translation to multiplication by the function ex{p). In 
view of the IV^- invar iance of Gq, the transformed function Q equals 

Q{p;z) = {:F,Q{z)r-^){p) = 5^Q(A,K;^)eA(p). (5.2) 

With the exponential estimate ()4.10j) of the function Q in mind, we infer that 
the function Q is well defined. It is also meromorphic in z with simple poles, 
which lie in a{Ho). Both the functions Q{p; z) and A{p) are real-analytic 
with respect to pi,p2 G M. 

Our goal now is to find the /^^-transformation of the Green function given 
by Krein's formula ()4.6|) . To this aim we denote 

RAiz) 

Roiz) 



= £ r 
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After a straightforward computation we arrive at the formula 

G{p; k, I, n; k', I', n; z) = 5kk' Snn' \ (5.3) 

e{l, n) — z 

-[Q{p;z) + A{p)] ^(y;;^ ,(^.,.)_, Xn(^3)XnK%), 
where So{p', k, I) is the /^^-transformed delta-function in M^, 

6,{p- k, I) = iV-i/2 ^ exp (2mm^^] ^*(7,Pi + m,l); (5.4) 

m=—oo ^ ' 

we have used here the relation 

Sx+^{p; k, I) = exp(7ri^ 7 A A) ex{p) S^{p; k, I) , (5.5) 

which follows from the VF,,-invariance of the free resolvent kernel and the 
fact that 6^{p;kJ) = {\B\{21 + 1) - 2)(£^Gg^(-, 7; ;z))(p; fc, /). Recall that 
although the Landau-Zak transformation £^ acts on L^(M^) £^(N), it can 
be viewed as two-dimensional since it does not affect the transverse modes. 

To perform the spectral analysis we need to know the behavior of the 
function Q{p] z) for fixed p and real z. It is convenient to treat this problem 
separately in each of the (infinitely many) intervals corresponding to gaps in 
the free Hamiltonian spectrum. To this end, we denote the points of (j{Hq) 
arranged in the ascending order by e^, i = 0, 1, ... . The function Q{p,-) 
diverges at a chosen point Si if and only if there exists at least one pair 
of integer numbers {l,n) G J{ei) such that Xn(^3) P belongs to 
Ui = {p & T"^ : 5o(p; - J) ^ 0}, where the expression 5o(p; ■, I) stands for an 
A^-dimensional vector. By ^ the residues in Krein's formula ()4.6|) are given 
by dQ/dz = Tl^T^; the transformed form of this relation reads 

Q/O 00 00 ^ iV— 1 

^fe = E E led n)-z\^ ^ 1^°^^' ^' (^-^^ 

/=0 n=l I V ) J I 

For notational convenience, we also put e_i = —00 and U-i = T^. It is clear 
from ()5.6|) that the function Q{p,-) is monotonously increasing in each of the 
intervals (ej-i, £j). 
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The asymptotic behaviour of Q{p, z) for large negative z is governed by 
the function Qol^s! z), while the contribution of the rest of the series in ex- 
pression ()5.2p for Q is bounded, its convergence being ensured by the estimate 
()4.1U|) . By [5J the divergent term has the expansion 

Qo{kz- z) = -^^^ + 0{l) as 2:^-00. (5.7) 

Now we are ready to analyze the spectrum of the fibre operator Ha (p) ■ The 
first thing we would like to know is whether the Landau levels of the free 
system stay in the spectrum a{HA{p))- By examining the residues of the 
Green function ()5.3|) at 2 = e a^Ho), = 0, 1, . . . , it is straightforward 
find the multiplicity of Ei in a{HA{p)) and the corresponding eigenspaces. 

Lemma 5.1 Assume that the flux r] equals an integer number N and the 
lattice r is 'monoatomic'. Fix a point Si from cr(i/o) o-nd p G T^, then one 
of the following three situations occurs: 

(i) if there is at least one pair {l,n) G J(ej) satisfying 6o{p;-,l) 7^ and 
Xnii^s) 7^ 0, the multiplicity d{p;ei) equals N\J{ei)\ — 1 and the eigenspace is 
the orthogonal complement of the vector uj{p; 6i) in the space fl{p]ei), where 



\ •^"^ / k,l,n 

(/o,no)eJ(£i) 

n{p;e,) = C^®( (Q®Xn)), (5.8) 

(l,n)eJ{ei) 

(ii) if Xn{i^3)So{p', -J) is a zero vector for all indices {l,n) G J{ei) and 
Q{p;ei) + A{p) 7^ 0, then d{p;ei) = N\J{ei)\ and Q{p]ei) is the eigenspace, 

(Hi) if Xn{i^3)So{p] ',^) is a zero vector for all indices {l,n) G J(£j) and 
Q{p; Ei) + A{p) = 0, then d{p; ei) = N\J{ei) \ + 1. In this case the eigenspace 
is the linear hull ofVLip^Si) and the vector l^°i£:Ml^A!lsl\ ^ where we put 

V '^^ / k,l,n 

0/0 := for elements with {l,n) G J{ei). 

This answers the question what happens with the free Hamiltonian spec- 
trum under influence of the perturbation. However, the spectrum a{HA{p)) 
contains not only the original modified Landau levels but also additional 
points due to the presence of the point potential. From the properties of 
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Q{p, ■) it is obvious that there is exactly one solution Ei^{p) of the implicit 
equation 

Qip;E) + Aip) = (5.9) 

in any interval (£jj,_^(p), £j^(p)), where {£ik(p))T=-i ^ subsequence of all points 
from {ei)°Z_i which the function Q{p; ■) diverges. If E does not belong 
to cr(iJo); the said solution is a nondegenarate eigenvalue of Ha{p) with the 
unnormalized eigenvector (^ ^g_°n)-Elp) Xn(%) j • This vector is non-trivial, 
which follows from the inequality derived below and applied to the 

derivative dQ/dz given in the monoatomic case by ()5.6|) . 

Let us summarize the effect of the periodic point potential on the spec- 
trum a{HA{p)) at a fixed point p G T^. The modified Landau level Ei in the 
spectrum of the free Hamiltonian Hq has the multiplicity equal to J(£:j)|. 
In the perturbed spectrum a{HA{p)) the generic situation is the case (i) of 
the above lemma, when the function Q{p; ■) diverges at Ei. Then an eigen- 
value splits off this level moving down with the increasing A{p) towards the 
neighbouring lower modified Landau levels. In particular, for N\J{Ei)\ = 1 
the perturbation removes in this way the level Ei from the spectrum a{HA{p)) 
entirely. On the other hand, if Q{p; ■) does not diverge at Ei, the multiphc- 
ity remains the same or it can be enlarged by one; the latter happens when 
an eigenvalue coming from a higher modified Landau level reaches Ei for a 
particular A{p). These two situations correspond, of course, to cases (ii) and 
(iii) of Lemma 15.11 respectively. 

After we have found the spectrum of the fibre operator for a fixed quasi- 
momentum, we can proceed to analysis of the spectrum of the full operator 

Ha^Ha= / HA{p)dp. (5.10) 

Since the functions appearing in the equation ()5.9|) are real-analytic the same 
is true for its solution. Fix an interval (ej_i,ej) such that there are two 
pairs of indices (/i,ni) G J^Ei^i) and (/2,'"-2) ^ J{.^i) with Xn^_{.i^z) 7^ and 
Xnii'^^) 7^ 0) then by the above discussion there is a real-analytic function 
Ei{-) defined on a set Ui^ fl Ui^ with the range in {Ei^i,Ei). Moreover, A{-) 
is bounded as a continuous function on a compact set. Combining this ob- 
servation with the asymptotics (|5.7|) for i = 0, we see that also the range of 
Ei{-) is a bounded interval. Since 5q{-] k, I) is an analytic function and Ui ^ 
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for all / G NU {0}, the domains Ui are dense open sets of full measure, so any 
intersection [/^^ fl Ui^ is also an open set of full measure. Hence the function 
Ei{-) extends by continuity to the entire torus and its range lies in the 
interval having a finite lower bound for i = 0. 

A modification is needed in case of an 'orphan' modified Landau level, 
i.e. a point ei> for which there is no pair of indices {l,n) G J{ei>) satisfying 
Xnif^s) 7^ 0. It is obvious that this cannot be a pole of Q{p; ■) for any p, and 
consequently, the implicit equation ()5.9|) may have no solution in one or both 
of the intervals and {ei',ei'+i). Instead we have to consider in this 

case the joint interval amended with the common endpoint. Then there is 
a unique solution of (j5.9|) on the interval (ej/.i, ^jz+i), provided its 

endpoints belong to the 'regular' class considered above, and the dispersion 
function Ei>{p) is excluded from further consideration; the argument easily 
extends to the situation with two or more neighbouring 'orphan' points. Note 
that the number of 'non-orphan' levels is infinite. 

Eliminating the 'orphan' points e^' from (£j)^o obtain a subsequence 
(£ij.)^Q; we add conventionally e^i = — oo as its first term. For each interval 
ejj.) we have then a unique dispersion function Ei^{p) defined on an 
open set of full measure as a solution to ()5.9p and extended by continuity to 
the entire torus T^. 

Lemma 5.2 The function Ei^{-) defined above has following properties: 

(i) if ei^_^ < Ei^ij)) < Si^, then Ei^^^p) is the unique solution to the implicit 
equation AS.yp in {ei^_-^,eij, 

(a) if is a pole of the function Q{p; ■), then Ei^{p) < < Ei^^^{p), 

(Hi) if Si^ is not a pole of the function Q{p]-), then 

Q{p; EiJ + A{p) < ^ E,^{p) = < E,^^^{p) , 

Q{p] Ei^) + A{p) > ^ E,^ (p) < Ei^ = E,^^^ {p) , 

Q(p;£iJ + I(p) = ^ E,^{p)=Ei^ = E,^^^{p). 

Proof: The proof is similar to the one for point interactions on the plane (7] . 
By definition of the sequence {Sif,}kL-i there are pairs of indices (/i,ni) e 
J(ej^_J and {l2,n2) G J{Ei^) such that Xnil^s) ^ and Xnal^a) 7^ 0. 
The union Ui^ fl Ui^ has a full measure in T^, so the extension of Ei^[-) 
to the whole by continuity (using a sequence G Ui^ fl Ui^ tend- 

ing to p G T^) is well defined. Assume first that Ei^{p) does not coin- 
cide with the endpoints and ej^. of the interval. From the condition 
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niPn, Ei^iPn)) ■= Q{Pn\ -^i^ (pn)) + ^(Pn) = and the joiiit continuity of ■) 
in the a neighbourhood of the point (p, Ei^{p)) we infer that Ei^{p)) = 
for any p E T"^; then the claim (i) follows from the monotonicity of the 
function ■). 

The point Ei^, (p) defined as above cannot be a pole of the function fi{p; ■). 
To see this, consider the function z z) := z){z — £:j^_J(2; — ej^.), 

which is analytic in an interval (^j^.^^ — p, + p) with a small enough p > 0. 
Using the continuity again we find that l3{p; Ei^{p)) = for any p & T"^, thus 
a pole at Ei^{j)) is excluded. This proves the claim (ii). 

To prove the last statement of the lemma, we need two auxiliary results: 

Eik (p) = Si^ ^ lim a{p;z) <0, 
Ei.ip) = ei,_, lim a{p; z)>0. (5.11) 

Let us check the first relation. Assume that Ei^{p) = but the limit is 
strictly positive, then there are Eq G {ei^^^.Si^) with fi{p; Eq) > and Pn 
such that fi{pn; Eq) > and Eq < Ei^{pn). However, the monotonicity of 
/J'ip;-) leads then to a contradiction, < fi{pn; Eq) < fi^pa; Eif,{pn)) = 0. 
Assume now that Ei^ is not a pole of fi{p; ■) and fi{p;ei^) < 0. By ()5.1ip 
we would have < Ei^^^{p). There is an integer number no such that 
jj,{pn',Si^) < holds for all n > no, hence Si^ < Ei^{pn) for all n > no, and 
consequently, Si^ < Ei^{p). The case of the opposite inequality nip^Si^,) > 
is treated in a similar way. 

Finally, let us consider the last case when Si^. is not a pole of p(p; ■) and 
p(p; EiJ = 0. We have to exclude both the strict inequalities Ei^{p) < and 

< Ei^^^{p). Assume, for instance, that the first one of them holds. From 
(jS.lip and the proof of the claim (i) we know that lim2^£;.^(p) p(p; 2;) > 0. 
There exists Ei satisfying (p) < Ei < Ei^, and we arrive at a contradiction, 
< fi{p; E,^{p)) <p(p;^i) <p(p;£^J =0. 

To finish the proof, we must check that the definition of Ei^, (p) is indepen- 
dent of the choice of the approximating sequence {pn)'^=i- Consider another 
sequence {p'^)'^^^ converging to p and denote the limit of Ei^^^p'^) by E[^{p). 
Assume Ei^^ij)) < E[^{p), then by the claim (i) it is necessary that at least one 
of these points coincides with one of the endpoints of the interval {Ei^_-^,Ei^). 
Using ()5.11|) we arrive at the relations 

lim u(»: > , lim u(n: z) < Q . 
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Choosing points Ei, E2 such that Ei^{p) < Ei < E2 < E[^{p), we get a 
contradiction, < Ei^{p)) < fi{p; Ei) < fi{p; E2) < n{p; E[^{p)) < 0. | 

Combining the above resuhs with with the direct-integral decomposition 
()5.10|) we arrive finally at the sought description of the spectrum a{HA)- 

Theorem 5.1 Suppose that the flux rj is integer, t] = N, and the elementary 
cell contains one point potential. Then a (Ha) consists of two parts: 

(i) The first one is the union of spectral hands denoted by Jk, = 0, 1, . . 
where Jk is the range of the function Ei^{-) over the torus T^, with Ei^{p) 
defined by the implicit equation i5.y\} . Each band Jk lies within one interval 
and two neighboring bands Jk and Jk+i have a common endpoint 
if and only if there exist pi andp2 from such that Q{pi; z) and Q{p2', z) 
do not have a pole at Si^ and 

g(pi;£,J + I(pi) > 0, 
Q{p2]ei,) + A{p2) < 0. 

There is at most one degenerate band corresponding to a constant Ei^{-). 
In particular, the degeneracy is excluded if the matrix A is diagonal, i.e. 
A = ale2(^r) for some a G M. The absolutely continuous spectrum of Ha is 
the union IJfclo '^i' with the exception of the possible degenerate band, 
(a) The point part of the spectrum consists, in addition to the mentioned 
degenerate band, of the modified Landau levels zq from (t{Hq) which persist 
under the perturbation. This concerns the whole (j{Hq) if N > 2 while for 
N = 1, the levels zq G ct^Hq) for which there is just one pair of indices 
{l,n) G J{zo) and Xnii^s) 7^ have to be removed. 

Proof: Most part follows from Lemma 15.11 and Lemma 15 .21 it remains us to 
check the claims about the degenerate band in (i). Suppose that there are 
two different degenerate bands {E} and {E'} with E, E' ^M. separated from 
the spectrum of Hq. Then we have Q{p] E) + A{j)) = for all p & T"^ and 
the same for E' which yields 

Q(7, K- E) = -A(7, k) = Q(7, K- E') (5.12) 

for all 7 from the set A. In particular, choosing 7 with a large modulus, we 
find different terms for the corresponding matrix elements of Q at different 
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energies E and E', because 



Q(7, /?; E) = C{E) e-¥H^|7|— (l + 0{\^r)) , (5.13) 



where u = min < 2, -[^(f )^ k which leads to a contradiction with ()5.12p . The 



expansion ()5.13p follows from the definition of Q by ()4.4j) in combination 
with the asymptotic behavior of the hypergeometric function, U{a,l;s) = 
s~"(l + - see 13.1.8]. Furthermore, consider a diagonal matrix 

A and suppose that there is a degenerate band {E}. Then the condition 
Q{p] E) + a = holds for any p G T^. In view of the relation ()5.2p it leads 
to the requirement QIj, k; E) = for all 7 7^ /? which again contradicts the 
known asymptotic behavior. | 

Thus we have obtained spectral bands between neighboring points of 
the unperturbed spectrum as in the planar case [7j. Needless to say, the 
bands are not the same because the dispersion functions are different and 
also the unperturbed spectrum is different: its points are sums of the 'two- 
dimensional' Landau levels and energies of transverse modes. In this sense 
the band structure in a layer is richer. 

The most important difference from the planar case is the possible exis- 
tence of a spectral gap containing the whole interval Ei) for some integer 
i, so that the free-spectrum gap is preserved by the perturbation. Such a sit- 
uation occurs, for example, if the positions of point potentials coincide with 
a node of each transverse mode corresponding to Ei, i.e. Xnii^s) = for all 
(/, n) G J{ei), and if at the same time Q{p; ei)+A{p) < holds for all p G T^. 
The last condition is satisfied, e.g., for a diagonal matrix A = al with the 
parameter a < — maXpgr2 \Q[p;ei)\ as it follows from monotonicity of the 
function Q{p] ■)■ Although it is not a generic situation it is not purely hypo- 
thetical. On the other hand, such a preserved gap cannot occur for the lowest 
interval, = 0, or in the situation when both endpoints contain the contribu- 
tion from the lowest transverse mode, i.e. (/, 1) G J(ei_i) and (/', 1) G J(£i). 
In case of a thin layer, it means that there is exactly one spectral band be- 
tween each two neighboring modified Landau levels below £:(0,2), i.e. below 
the threshold of the second transverse mode; it is obvious that the thinness 
here has to be understood in comparison to the characteristic length given 
by the magnetic field. 
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6 The case of integral flux and polyatomic 
crystal 

Consider next a polyatomic lattice A, i.e. suppose that the set K contains 
more than one point, Ri 7^ Rj for i, j = 1, . . . , \K\, i 7^ j. The flux 77 is again 
an integer number N . Compared to (j5.ip the Fourier transformation must 
be modified; : ^^(r) ^ L^{T^) ® f{K) acts now as 

(J^^0)(j9; /?) = ^ 0(A) exp [ni^ kAX- 2m{XaPi + XbP2 + NXaXb/2) ] , 
AeA 

(6.1) 

i.e. values of the transformed function JF^0 of the variable p E are no 
longer scalar but rather |i^r|- dimensional vectors. The Fourier transformed 
operator J^^D{{X, 1))JF~^ with A G A acts again as a multiphcation by e\{p). 
An argument similar to that of the monoatomic case leads to 

{Q{z) + A){p; K, k') = (J^rjiQiz) + A)J^-'){p; k, k) (6.2) 

= "^(Qiz) +A){X + R') exp [TTi^/t A A - 2m{XaPi + Xi,p2 + iVAaA6/2)] . 
aga 

The matrix elements Q{p:, R, R'; z) are real-analytic with respect to p G and 
meromorphic in z with simple poles which can be located only at the points 
of cr(ifo)- The matrix elements A{p] R, R') are real-analytic in p. Finally, the 
transformed Green function reads 

G{p] k,l,n] k', I', n; z) = 6kk'Sii'6nn' \ (6.3) 

e[l, n) — z 

- Y: [Q(p; z) + Aip)]-\R, R') ^'^ Xn{ns)Xnin',) . 

eiLn) — z e[l'n') — z 

k,k'&K \ ^ / \ ^ J 

Using the known relation dQ/dz = T**rz for the derivative, we get 
— (p; R, R'; ^) = X] I (I n)_ ^12 Xn{i^3)Xn{K3) "Y 5*(p; k, l)6^\p] k, I). 

^ 1=0 n=l l^*^ ' ^1 A;=0 

(6.4) 

The asymptotic behaviour of the diagonal elements of the matrix Q{p; z) as 
z —00 is similar to ()5.7|) . while the non-diagonal elements are bounded. 



Q{p;R,R';z) = -5^^,^ + Oil). (6.5) 
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Let us begin the spectral analysis with the fiber operator Ha{p) for a 
fixed p G T^. For each point Ei G (j{Hq) we define matrix G^Xp) ^is the 
residue of —Q{p; z) at z = Ei, i.e. 

G^Xp'^k,k')= ^ Xnif^3)Xnif^3)^S*^ip;k,l)6^'{p;kJ) , (6.6) 

(«,n)GJ(£i) fc=0 

and we denote its rank by r^. (j?) . We further define P^. (p) as the orthogonal 
projection onto keiGeXp) i'^{K) and operator Di;.{p) as 

D,Xp) = lim PeAp)iQ{P; + ^(p))|kerG.^(p) • (6.7) 

The multiplicity of Ei in the spectrum of the free Hamiltonian is equal to 
N\J{Ei)\ with the eigenspace VL^Ei) defined by (jS.Sj) . The second term of the 
Green function ()6.3p modifies the residue at z = Ei in two possible ways: 

(i) taking into account contribution from the indices (/,n), {l',n') G J{Ei) to 
the kernel of the resolvent, we get the orthogonal projection to the subspace 
spanned by the vectors (Xni('^3)^K(p; ^, ^i))fc,i with k = 0, . . . , N — 1 and 
{li^rii) G J{E.i). In this way the multiplicity of Ei is diminished by ^^.(p). 

(ii) On the other hand, for indices (/,ti) ^ J{Ei) we get a nonzero residue 
when the operator D^-{p) is not invertible. The corresponding eigenspace is 
orthogonal to VL{Ei) with the maximal possible dimension equal to \K\. 

Putting two terms of the Green function ()6.3|) together, we arrive at the 
following result: 

Lemma 6.1 Assume that the flux t] equals an integer number N and the el- 
ementary cell of the lattice T contains \K\ point interactions. Choose a point 
Ei G (j{Hq) and fix p G T^. Then the multiplicity ds-{p) of Ei in a{HA{p)) is 
equal to 

4,(p) = iV|J(e,)|-r,,(p), 
if Di;-[p) is invertible, while in the opposite case it satisfies the inequalities 

N\JiE.)\ - r^Xp) < dsM < N\J{e::)\ + \K\ - r^Xp) ■ 

Apart from the modified Landau levels of the free system, the spectrum 
of Ha{p) contains eigenvalues due to the presence of the point potentials. A 
necessary condition for G M \ (t{Hq) to be an eigenvalue is 

det[Q(p;E) + I(p)] = 0. (6.8) 
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> cJp(^K) (6.9) 



Notice that E might not be an eigenvalue of Ha{p) if the vectors 

\e[l,n) — E / k,i,n 

were not hnearly independent, in which case the second term in Green func- 
tion ()6.3p could vanish. This cannot happen, however, because dQ/dz is 
the Gram matrix for this |fr|-tuple of vectors and by one has for a fixed 
z ^ (j{Ho) the inequality 

dQjp] z) 
dz 

with some Cz > 0. Therefore ipft are linearly independent and eigenvectors 
corresponding to E are given by J^geK f^f^'^f^^ where the vectors belong 
to ker[Q(p; E) + A{p)]. 

A useful way to solve the implicit equation ()6.8|) is by examining the 
eigenvalues jJi^^\p] z), j = 1, . . . , \K\ of the matrix Q{p; z) + A{p) with the 
numbering which takes their multiplicity into account. Apparently, a solu- 
tion of equation fx^-'^p; E) = for some j solves also the original equation. 
Properties of the functions /i are described in following lemma: 

Lemma 6.2 Suppose that the flux r] is an integer number N and the ele- 
mentary cell of r contains \K\ point potentials, and fix p G T^. Then the 
eigenvalues f^^^\p; z) of the matrix Q{p; z) + A{p) are monotonously increas- 
ing functions of z in each interval {si^i^Si), z = 0, 1, . . . . When z approaches 
a modified Landau level Si, exactly rg.{p) functions of the family n^^\p; ■), 
j = 1, . . . , li^l diverge, while the remaining ones have finite limits. When z 
tends to —oo, all the \K\ functions fJ,^^\p; ■) diverge to — oo. 

Proof: Since the matrix Q{p;z2) — Q{p;zi) with < Zi < Z2 < Ei is 
by ()6.9|) positive definite, the Lidskii's theorem (TUl Thm. II. 6. 10] yields the 
inequality iJ.^^\p; zi) < iJ.^^\p; Z2) for each j, hence all the functions iJ,^^\p; ■) 
increase monotonously between every two neighboring points of (j{Hq). To 
prove the second claim, take \K\ eigenvalues P^^\p;z), j = 1,...,\K\, of 
matrix {si — z)[Q{p;z) + By Rellich's theorem ^01 Thm. II. 6. 8] the 

functions P^^^ are continuously differentiable in some neighborhood of Ei, 
an for z = the matrix under considerations coincides with Gg.. Finally, 
using Gershgorin's circles - see [HI Thm. XIV. 5. 5] - in combination with the 
asymptotic formula (j6.5j) for the matrix Q{p; z) we get the last claim. | 
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The lemma implies that there are at most \K\ eigenvalues in each interval 
Ei), i = 0,1, .. . which is after all clear from general principles - cf. 
Sec. 8.3]. To get a more specific information about the actual number of 
eigenvalues, we have to look closely at the number ^^^(p). 

Remark 6.1 Since the matrix G^Xp) is obviously the Gram matrix of the 
system of J(£:j)|-dimensional vectors 

[Xns{i^3)^^{p;k,ls)) k = 0,...,N -1, {ls,n,) E J{ei) , 

\ / k,s 

the rank r^.{p) cannot exceed min(|i^|, J(£:j)|). There are three situations 
when the maximum value specified here cannot be achieved whatever point 
p is considered. First, the dimension of the vectors is in fact smaller if 
Xusii^s) = holds for some index Ug and all k E K. In general, the said 
dimension is equal to N\J[ei)\, where J{ei) is the set J{ei) from which the 
pairs {ls,ns) with the described property were deleted. Second, the number 
of the vectors is li^l, but when a lattice point /? is such that Xnsl'^s) = holds 
for all {ls,ns) G J{si), the corresponding vector is zero for any p. Excluding 
such elements from the set K we obtain a subset denoted as K{ei). 

Finally, we must also examine carefully the situation when several points 
, j = 1, . . . ,q, from the elementary cell differ in the third coordinate only, 
i.e. they are arranged vertically in the layer. If the level Ei is not degenerate, 
by neglecting all but one of the q vectors in question, we do not change the 
rank of the Gram matrix, provided the one which we keep is not a zero vector. 
If a degenerate level is admitted, then the number of linearly independent 
vectors among these q vectors is less or equal to the rank of the following 
matrix 

(xnsii^i)) j = 0,...,q, {ls,ns) e J (Si) . 

Then we eliminate the remaining ones of the q vectors from K{ei) and obtain 
in this way another subset of K which we denote as K{ei). The maximal 
possible rank r^. (p) of the matrix . (p) is in this way equal to 

mmi\Kie,)\,N\Jie,)\); 

whether the maximum is achieved at a point p or not now depends on the 
functions k, I). 

Such a dependence of r^. {p) on the parameters of the problem makes the 
general spectral analysis cumbersome. In what follows we will thus restrict 
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our attention to the generic situation only and we impose additional restric- 
tions on the Hamiltonian H^- First of all, we assume that {n/df and \B\ 
are not rationally related, so that J{ei) = 1 holds for all modified Landau 
levels Ei. Furthermore, we define the sets 

:= {peT^: r,^(p) = w = min(iV, |^|)} , 
U'J^ ■= {peT^ : D,^ is invertible} 

and 

f/; n f/;; , (6.10) 

where \K\ represents the number of points k & K after we have excluded 
q — 1 points from every g-tuple in which the first two coordinates coincide. 
The set K does not depend on Si any more. 

In the rest of the paper, we consider only operators Ha such that the 
set nj>o nonempty, which is a generic situation. One can check easily 
that by an arbitrarily small shift of the points k in the elementary cell Qa 
we can achieve that ni>o ^'^'-^ ^ similar way, by an arbitrarily 

small perturbation of the diagonal elements of the matrix A we can always 
satisfy the condition ni>o ^e' ^- Recall that the possible positions k of 
point potentials in the elementary cell are dense in Q^^\ while all possible 
values of A{k,, k) span R'^', and consequently, the set 

\ (^1, ...,K\K\, A{ki, Ki),..., i^iK])) e gf X rI^I : fl Us,, + \ 

I i>0 J 

has a full measure in Qjf' x R'''^'. 

For \K\ = Tmax and p G f/ei„i fl t/g-, the number of the eigenvalues in the 
free Hamiltonian spectral gap equals \K\. Now we employ the assumption 
C\i>oU£i 7^ ^- since S^ip^kJ) and Q{p;z) with z ^ cr{Ho) are analytic as 
functions of p, all sets the [/^. and U'J., i > 0, have full measures, and 

the same is true for the intersection f/e._j fl f/^-. The functions e\^\-) are 
bounded, which follows for i = from the asymptotic formula ()6.5|) and 
the boundedness of the matrix A while otherwise the claim is valid trivially. 
Thus we may extend the \K\ dispersion functions e\'^\-) by continuity to the 
entire torus T^. We arrive at the following result. 
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Lemma 6.3 Assume that \K\ = rmax (^nd the Hamiltonian Ha satisfies 
C\i>oUe, 7^ 0- Consider the extended functions e\^\p), j = 1, . . . , \K\, de- 
fined on in the described way. Then the inequalities Si-i < e\^\p) < Si 
for some j, 1 < j < l-ft'l, implies that E^''\p) is a solution of the implicit 
equation \d. t^) . 

For li^l > rjnax and p G f^^.i n 11^^, the number of eigenvalues in the 
interval is not necessarily equal to li^l. Consider first p G nj>o^£i- 

In this case, \K\ — rmax functions fj,^^\p; z) do not diverge when z approaches 
an endpoint of the interval, instead they 'meet' there one of functions from 
the neighboring interval. Then it is natural to unify these functions coming 
from neighboring intervals obtaining new functions i'\^\p] ), i = 0, 1, . . . , j = 
1, • • • 5''^max5 which are defined on the enlarged intervals {Si^s^Si), where 1 < 
s < l-^l/^max if the fraction is integer, otherwise 1 < s < [l-ft'l/rmax] + 1- We 
put Ei-s = — oo and t/e,.^ = whenever i — s < 0. On each subinterval, the 
corresponding families of functions z/(p; z) and fi{p; z) coincide. The function 
v\''\p;z) tends to ±oo as z approaches Si or ei-s, respectively. Thus there 
exists exactly one solution to the equation v'f \p]E) = which we denote 
as The same conclusion can be made under a weaker assumption, 

p e riLo ^-^i-r '^^^^ again a full measure and functions e\^\-) are 

bounded, so we can extend the dispersion functions by continuity to the 
entire torus T^. 

After having analyzed the fiber operator spectrum, let us run p run 
through the torus to get the spectrum of the Hamiltonian Ha- 

Theorem 6.1 Suppose that the flux rj is an integer number N and the ele- 
mentary cell of the lattice T contains \K\ points. In addition, let {ir/dy and 
\B\ be irrationally related and let the Hamiltonian Ha satisfy the condition 
C\i>o 7^ ^' "^here the sets 1/^^ are defined by i6.1(J\} . Then the spectrum of 
Ha consists of two parts, namely 

(i) spectral bands with j = 1, . . . , rmax and i = 0, 1, ... , where rmax = 
min(A^, \K\) and K is the maximal subset of K such that no pair of points 
from K' coincide in the first two coordinates. Each band Ij is given as the 
range of the extended function e\^\p), p G T^, which is defined by the implicit 
equation \6. it lies within the interval [si-s^Sj], where 1 < s < li^l/rmax if 
the fraction is integer and 1 < s < [l-ft'l/rmax] + 1 otherwise. The absolutely 
continuous spectrum of Ha is the union IJi^o Uj=r H except possible bands 
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degenerate to a point, 

(a) modified Landau levels from the spectrum of Hq. If N < rmax? some 
points of ct^Hq) may be absent from the spectrum of Ha. 

The theorem says nothing about possible common endpoint of two neigh- 
boring bands, which coul be compared to the analogous part of Theorem 15.11 
in the monoatomic case. However, the term 'neighboring bands' does not 
have much sense here; the bands // with the same index i may overlap and 
for \K\ > Tmax also bands with different indices i may overlap. 

Apart from the difference between the modified Landau levels here from 
the unperturbed spectrum in the planar case, the structure of the two spectra 
is similar with one difference: the number of bands // which neighbor with the 
same point Ei from above equals min(A^, \K\), while in the planar case it 
is min(A^, \K\) instead. The reason is clear: the magnetic field perpendicular 
to the layer does not 'distinguish' two points placed one on the top of the 
other, and such a situation can never occur in the planar case. 

In the previous chapter we have found that the spectrum can contain a gap 
covering the whole interval {ei-i,ei) for some i. In the polyatomic case it is 
obviously possible for l-ft'l > 2rmax5 while the above discussion excludes such 
a situation otherwise. It might thus seem that for \K\ = rmax = 1 we get a 
contradiction. However, the discrepancy comes from the stronger restriction 
we have imposed upon the operator Ha] the condition nj>o 7^ ^ does not 
allow Xnl'^s) = to hold for all (/,n) G J{e). 

7 The case of a rational flux 

The general case when the flux is a rational number, V = jj^ can be in 
some sense reduced to the previous analysis. We can pass to an integral flux 
7]' = N hj enlarging the elementary cell. The new lattice A' is generated by 
vectors a and Mb and the new set K' is given by + {0, b, . . . , (M — l)b}. 

However, the result obtained by this simple trick is not fully correct be- 
cause it does not take the M-fold degeneracy into account. Recall that the 
Hamiltonian Ha commutes with all magnetic translations from Wrj. Hence 
any of its eigenspaces can be written as a direct sum of spaces of the irre- 
ducible representations of Wrj, which are M- dimensional. Therefore addi- 
tional modifications are needed here. We define the Fourier transformation 
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jr^ : £2(r) ^ L2(r2) ® ® ® ^^(is:) by the prescription 
{J='r,(j)){p;j,m,K) = 0(Aaa+ (A6M + m)6 + /«) 

X exp K A (Aaa + Xj^Mh) 



(7.1) 



X exp 



m + 2j 
M 



The transformed quantity J^r)D[-)T^ ^ is then a direct integral of multiples 
of irreducible representations, 

The representation Ad(-;p) acting on C-^ ® C-^ is given by 

Ad((a,l);p) = A((a,l);p)®/cM, 
A,((M);p) = 5 0A'(p), 

where 5" and A'(p) are operators on which act at the basis vectors e^, 
J = 0, . . . , M — 1 , in following way, 

Sej = e^ei , A'(p)ej = exp(-27rip25i,M-i)ej®i , 

with © and representing the sum and difference in {0, . . . , M — 1} modulo 
M. This representation is equivalent to the representation A(-;p) ® Ic^i . 
Applying now the Fourier transformation to the operator Q{z) + A we 

get 



{Q{z) + A) (p; j', m, j', m', = 5jji exp vri^ k A (AqO + XhMh) 

iV 



X exp 



n • / X X X /x m + 2j\ 

-27rz I AaPi + XhP2 + y ''^a l^'^fe ) 

X (Q(^) + A){\ad+ {XhM + m)b + K,m'b + k') . 
Finally, the transformed Green function reads 

G{p;j,k,l,n;j',k',l',n';z) = SjfSkk'Sii'd„ 



1 



e(l, n) — z 



(7.2) 



(7.3) 



27 



M-1 
K,K.'£K m,m'=0 

^ eil,n)-z eil',n')-z ^"("3)Xn^(^3) , 

where 

~ °° f P2 + k\ 

SmAP'J^kJ) = N^'^^'^ ^ exp ( 2TTir j i/j^imb + K;pi + r]j + rj) . 

r=—oo ^ ' 

(7.4) 

We have employed a relation similar to ()5.5p with A G A' and an additional 
parameter j, and the M-dimensional representation A(-,p) instead of the 
character xi'iP)- 

Next, we proceed to the spectral analysis in the same way as we did in 
the previous section. The implicit equation defining the dispersion function 
{p) has the form 

dei[Q{p-j-E)+A{p-j)] = Q (7.5) 

with the dimension of the matrix being equal to |i^|M. Due to the equiva- 
lence of the two representations which we have mentioned above, the matrices 
Q{p;j; z) + A{j)]i) for different j are unitarily equivalent. This is the source 
of the M-fold degeneracy of the eigenvalues of Ha{p), and thus also of the 
spectral bands of Ha- 

To simplify the description of the spectrum a{HA), we consider in the 
following theorem a lattice with \K\ = 1. 

Theorem 7.1 Assume that the flux rj is a rational number N/M and the 
elementary cell contains one point potential. In addition, let (vr/d)^ and 
\B\ be irrationally related, and let the Hamiltonian Ha satisfy the condition 
f]j>Q r\f=^ 7^ where the sets 11^. are defined in analogy with the expres- 
sion id.lU^) . Then the spectrum a{HA) consists of two parts: 
(i) The first one is the union of spectral bands with r = 1, . . . ,rmax o-nd 
i = 0, 1, . . where rmax = min(M, A^). The band II is given as a range of the 
function p t-^ Ef'^\p) with p G T^, defined by the implicit equation \7.5^ for 
some j . Each band II is M -times degenerated and it lies within the interval 
[si^s.^^i], where l<s<MifN = l and 1 < s < [M/N] + 1 otherwise. The 
absolutely continuous spectrum is the union of all bands II except possible 
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bands degenerated to a point. 

(a) The second part of the spectrum contains modified Landau levels from 
a{Ho). In the case N < M some points of a{HQ) may he absent from the 
spectrum a {Ha) ■ 

8 Survey of the results 

We have analyzed the spectrum of a Dirichlet layer with a periodic array of 
point perturbations in presence of homogeneous magnetic field. The generic 
picture we have obtained for a rational flux, r] = N/M, has some well-known 
features analogous to [7j: in the case of a single potential in the elementary 
cell, there are min(M, A^) spectral bands which split off each modified Landau 
level Ei, each band is M-times degenerated, and its location is not necessarily 
restricted to the gaps adjacent to the Landau level in question; more precisely, 
if M > it may spread below ^i-i- In the case of n > 1 perturbations in the 
elementary cell, the number of the spectral bands changes to min(nM, A^), 
while the M-fold degeneracy remains the same. 

Apart of the magnetic field and the lattice spacing, the system has another 
parameter, namely the layer width d. Its first and most visible effect on the 
spectrum is that the Landau levels in the unperturbed spectrum are combined 
with the energies of transverse modes (that is what we mean by modified 
Landau levels), and thus they are described by a pair of quantum numbers. 
Due to this fact the results are similar to those for the two-dimensional system 
of j7j only in the generic situation when {TT/d)"^ /\B\ is an irrational number 
and no point potential is placed at a node of a transverse mode. 

If these additional conditions are not satisfied, one has to examine each 
Landau level separately as we did in Remark l6.ll For example, if (vr/d)^ and 
\B\ are rationally related, some Landau levels have an extra degeneracy and 
the number of bands may increase, because one must consider appropriate 
multiple of the integer A^. On the other hand, if we allow the sites of point 
potentials coincide with a node of a transverse mode, the number of bands 
may decrease because unperturbed levels do not "feel" the interaction. 

Acknowledgement 

The authors thank V.A. Geyler for useful comments. The research was par- 
tially supported by GAAS grant A1048101. 



29 



References 

[1] M. S. Abramowitz, 1. A. Stegun, eds.: Handbook of Mathematical Func- 
tions, Dover, New York 1965. 

[2] S. Albeverio, F. Gesztesy, R. H0egh-Krohn, H. Holden: Solvable Models 
in Quantum Mechanics, Springer, Heidelberg 1988. 

[3] V. V. Dodonov, I. A. Malkin, V. I. Man'ko: Phys. Lett. A51, 133 (1975). 

[4] P. Exner, R. Gawlista, R Seba, M. Tater: Ann. Phys. 252, 133 (1996). 
[5] P. Exner, K. Nemcova: J. Math. Phys. 43, 1152 (2002). 
[6] F. R. Gantmacher: The Theory of Matrices, Nauka, Moscow 1966. 
[7] V. A. Geyler: St. Petersburg Math. J. 3, 489 (1992). 

[8] V. A. Geyler, V. A. Margulis, I. I. Chuchaev: Sibirsk. Mat. Zh. 36, 828 

(1995). 

(English translation) Siberian Math. J. 36, 714 (1995). 
[9] V. A. Krein, G. H. Langer: Funktsional. Anal, i Prilozhen. 5, 59 (1971). 

[10] T. Kato: Perturbation theory for linear operators. Springer, Berlin 1976. 

[11] W. Opechowski, W. G. Tarn: Physica 42, 529 (1969). 
[12] A. Posilicano: J. Func. Anal. 183, 109 (2001). 

[13] J. Weidmann: Linear Operators in Hilbert Space, Springer, New York 
1980. 

[14] J. Zak: Phys. Rev. 134 A1602 (1964). 



30 



